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(57) Abstract 

The invention relates to an improved adaptive spectral estimator for estimating the spectral components in a signal (205) containing 
both an information signal, such as speech, and noise. A method and system provide for generating noise estimates (245) and then 
only updating the noise estimates during pauses in an information signal, when speech or other information is not detected, rather than 
continuously updating the noise estimates. A noise estimate is calculated for each frequency band and provides for the inclusion of a 
variable mathematical factor that can be set by the user to produce the best sound quality. 
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1 A METHOD AND SYSTEM FOR UPDATING NOISE ESTIMATES DURING 

2 PAUSES IN AN INFORMATION SIGNAL 

3 

4 FIELD OF THE INVENTION 

5 

6 This invention relates to a method and system for improving Adaptive Speech 

7 Filter (ASF) estimates of the noise component of complex signals that contain both the 

8 information signal and noise. The present invention generates noise estimates that are 

9 updated only during pauses of the information signal. This produces an increase in 

10 processing speed and a decrease in system memory. The methods of the present 

11 invention are particularly suited to implementation on inexpensive digital signal 

12 processors. 
13 

14 BACKGROUND OF THE INVENTION 

15 

16 The spectral components of an information signal are used in a number of signal 

17 processing systems including channel vocoders for communication of speech, speech 

18 recognition systems and signal enhancement filters. Since the inputs to these systems 

19 are often contaminated by noise there has been a great deal of interest in noise reduction 

20 techniques and consequently noise estimation techniques. The effect of uncorrelated 

21 noise is to add a random component to the power in each frequency band, and the 

22 subject of accurately assessing the noise content is crucial to achieve the desired end 

23 result, which is the elimination of noise from the complex signal. 
24 

25 Noise-free spectral components are required for optimum operation of channel 

26 vocoders. In a vocoder the input signal is filtered into a number of different frequency 

27 bands and the signal from each band is rectified (squared) and smoothed (low pass 

28 filtered). The smoothing process tends to reduce the variance of the noise. Such 

29 methods are disclosed in U.S. Patent No. 3,431,355 to Rothauser et al and U.S. Patent 

30 No. 3,431,355 to Schroeder. An alternative approach is disclosed in U.S. Patent No. 

31 3,855.423 to Brendzel et al. In this approach the level of the noise in each band is 

32 estimated from successive minima of the energy in that band and the level of the signal 
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2 

1 is estimated from successive maxima. In U.S. Patent No. 4,000,369 to Paul et al, the 

2 noise levels are estimated in a similar fashion and subtracted from the input signals to 

3 obtain a better estimate of the speech signal in each band. This method reduces the 

4 mean value of the noise. 
5 

6 Another application of spectral processing is for speech filtering. Weiss et al., in 

7 "Processing Speech Signals to Attenuate Interference", presented at the IEEE Symp. 

8 Speech Recognition. April 1974, disclose a spectral shaping technique. This technique 

9 uses frequency domain processing and describes two approaches - amplitude modulation 

10 (which is equivalent to gain control) and amplitude clipping (which is equivalent to a 

11 technique called spectral subtraction). Neither the noise estimate nor the speech 

12 estimate is updated so this filter is not adaptive. An output time waveform is obtained 

13 by recombining the spectral estimates with the original phases. 
14 

15 An adaptive speech filter is disclosed in U.S. Patent No 4,185,168 to Graupe and 

16 Causey, which is included by reference herein. Graupe and Causey describe a method 

17 for the adaptive filtering of a noisy speech signal based on the assumption that the noise 

18 has relatively stationary statistics compared to the speech signal. 
19 

20 in Graupe and Causey's method the input signal is divided into a set of signals 

21 limited to different frequency bands. The signal to noise ratio for each signal is then 

22 estimated in accordance with the time- wise variations of it's absolute value. The gain of 

23 each signal is then controlled according to an estimate of the signal to noise ratio (the 

24 gain typically being close to unity for high signal to noise ratio and less than unity for 

25 low signal to noise ratio). 
26 

27 Graupe and Causey describe a particular method for estimating the noise power 

28 from successive minima in the signals, and describe several methods for determining the 

29 gain as a function of the estimated noise and signal powers. This is an alternative to the 

30 method described earlier in U.S. Patent No. 4,025,721 to Graupe and Causey, which 

31 detects the pauses between utterances in the input speech signal and updates estimates of 

32 the noise parameters during these pauses. In U.S. Patent No. 4.025,721, Graupe and 
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Causey describe the use of Wiener and Kalman filters to reduce the noise. These filters 
can be implemented in the time domain or the frequency domain. 

Boll, in "Suppression of Acoustic Noise in Speech using Spectral Subtraction", 
IEEE Transactions on Acoustics, Speech and Signal Processing. Vol. ASSP-27, No. 2, 
April, 1979, describes a computationally more efficient way of doing spectral 
subtraction. In the spectral subtraction technique, used by Paul, Weiss and Boll, a 
constant or slowly varying estimate of the noise spectrum is subtracted. However, 
successive measurements of the noise power in each frequency bin vary rapidly and only 
the mean level of the noise is reduced by spectral subtraction. The residual noise will 
depend upon the variance of the noise power. This is true also of Weiss's spectral 
shaping technique where the spectral gains are constant. In Graupe's method the gain 
applied to each bin is continuously varied so that both the variance and the mean level of 
the noise can be reduced. 

There are many schemes for determining the spectral gains. One scheme is 
described by Ephraim and Malah in "Speech enhancement using a minimum mean- 
square error short-time spectral amplitude estimator", IEEE Transactions on Acoustics, 
Speech and Signal Processing, Vol. ASSP-32, No. 6, Dec. 1984. This describes a 
technique for obtaining two estimates of the signal to noise ratio - one from the input 
signal and one from the output signal. It does not update the estimate of the noise level. 
The gain is a complicated mathematical function of these two estimates, so this method 
is not suitable for direct implementation on a digital processor. 

In U.S. Patent No. 5,012,519 to Aldersburg et al the gain estimation technique of 
Ephraim and Malah is combined with the noise parameter estimation method disclosed 
in U.S. Patent No. 4,025,721 to Graupe and Causey to provide a fully adaptive system. 
The mathematical function of Ephraim and Malah is replaced with a two-dimensional 
lookup table to determine the gains. However, since the estimates of the signal to noise 
ratio can vary over a very large range, this table requires a large amount of expensive 
processor time and memory. Aldersburg et al use a separate voice detection system on 
the input signal which requires significant additional processing time. 
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1 There is thus an unmet need in the art to be able to utilize an efficient adaptive 

2 signal processing technique for the accurate and fast identification of noise. Processing 

3 time and memory efficiency would be improved if the noise estimates were only done 

4 during pauses of the information signal, so that noise estimates are updated only when 

5 an information signal is not detected. The algorithm should be capable of being 

6 implemented on inexpensive digital signal processors. 

7 

8 SUMMARY OF THE INVENTION 

9 

10 it is an object of the present invention to be able to obtain and update noise 

11 estimates only during pauses of the information signal, thereby decreasing processing 

12 time and memory requirements. 
13 

14 Therefore, according to the present invention, a method and system provide for 

15 noise estimates to be updated only during pauses in an information signal, when speech 

16 or other information is not detected, rather than continuously updating the noise 

17 estimates. Waiting for pauses in the information signal before updating the noise 

18 estimates allows processing time and memory requirements to be decreased. It also 

19 allows adaptive speech filtering to be easily implemented on inexpensive digital signal 

20 processors. 
21 

22 According to the method of the present invention, after a set of input frequency 

23 components have been produced, the total power calculated for each input frequency 

24 component, the power of the information signal estimated, a modified gain of the 

25 information signal calculated, and the input frequency component multiplied by the 

26 modified gain to produce an estimate of the power of the frequency component, then an 

27 estimate of the noise power is updated only if a pause in the information signal has been 

28 detected. Detecting a pause in the information signal is accomplished by first 

29 determining whether the estimate of the power of the frequency component of the 

30 information signal exceeds a first predetermined threshold value at each frequency. If 

31 the estimate of the power does exceed the first predetermined threshold value, then a 

32 threshold value thrsholdCnt[f] is checked to determine if it exceeds a second 
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1 predetermined threshold value. If the threshold value does exceed the second 

2 predetermined threshold value, then a pause has been detected. If the threshold value 

3 does not exceed the second predetermined threshold value, then no pause has been 

4 detected. In this instance, the noise estimate is not updated and instead the threshold 

5 value thrsholdCnt[f] is incremented. 
6 

7 The foregoing method of the present invention is implemented by a system for 

8 estimating the noise power of frequency components of an information signal from an 

9 input signal containing both the information signal and noise. The system has means to 

10 produce input frequency components, one frequency component for each frequency 

11 band, a first calculating means for calculating the total power of each input frequency 

12 component, a second calculating means for calculating the modified gain of each 

13 frequency band, and a gain multiplying means for multiplying the input frequency 

14 component by the gain to produce an estimate of the power of the frequency component 

15 of the information signal. The system additionally has an estimating means that 

16 estimates the power of the information signal and updates the estimate of the noise 

17 power only during a pause detected in the information signal by the estimating means. 

18 The estimating means itself has an adder, a first comparison element coupled to the 

19 adder that receives the estimate of the power of the information signal and a first 

20 predetermined threshold value from the adder and determines whether the estimate of 

21 the power of the information signal exceeds the first predetermined threshold value, and 

22 a second comparison element coupled to the first comparison element that determines 

23 whether a threshold value exceeds a second predetermined threshold value if the 

24 estimate of the power of the information signal exceeds the first predetermined threshold 

25 value. The estimate of the noise power is updated if the threshold value exceeds the 

26 second predetermined threshold value since this condition is indicative of a pause 

27 detected in the information signal. If a pause is not detected then an increment element 

28 increments the threshold value. 
29 

30 
31 



I 
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! BRIEF DESCRIPTION OF THE DRAWINGS 

2 

3 The novel features believed characteristic of the invention are set forth in the 

4 appended claims. The invention itself, however, as well as a preferred mode of use, and 

5 further objects and advantages thereof, will best be understood by reference to the 

6 following detailed description of an illustrative embodiment when read in conjunction 

7 with the accompanying drawing(s), wherein: 
8 

9 Figure 1 is a block diagram of a prior art system; 
10 

{ ! Figure 2 is a block diagram of a typical system of the present invention; 
12 

13 Figure 3 is a block diagram of a sub-system for gain modification, according to 

14 the present invention; 
15 

16 Figure 4 is a block diagram of a sub-system for signal power estimation, 

17 according to the present invention; 
18 

19 Figure 5 is a block diagram of a sub-system for noise power estimation, 

20 according to the present invention; and 
21 

2 2 Figure 6 is a block diagram of a sub-system for an information signal detector, 

23 according to the present invention. 
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! DESCRIPTION OF THE INVENTION 

2 

3 The present invention describes a method for generating noise estimates which 

4 are only updated when an information signal, such as speech, is not detected. The noise 

5 estimate is calculated for each frequency band, and provides for the inclusion of a 

6 variable mathematical factor that can be set by the user in order to produce the best 

7 sound quality. The noise estimate method of the present invention allows the adaptive 

8 speech filter algorithm to perform better under all conditions. 
9 

10 The adaptive speech filtering of the present invention is a modified version of 

11 that described in U.S. Patent No. 4.185,168 to Graupe and Causey which describes a 

12 method for the adaptive filtering of a noisy speech signal. The method is based on the 

13 assumption that the noise has relatively stationary statistics compared to the speech 

14 signal. 
15 

16 The input to the filter is usually a digital signal obtained by passing an analog 

17 signal, containing noise and the information signal, through high- and low-pass filters 

18 and then sampling the resulting signal at a sample rate of at least 8 kHz. The high pass 

19 filter is designed to remove low frequency noise that might adversely affect the dynamic 

20 range of the filter. The turnover frequency of the high pass filter is less then fjow , 

21 where fjow is the lower limit of the speech band in Hertz. The low pass filter is an 

22 anti-aliasing filter, which has a turnover frequency of at least fjiigh , where fjiigh is 

23 the upper limit of the speech band in Hertz. The order of the low pass filter is 

24 determined by the sampling frequency and the need to prevent aliasing. 
25 

26 The output signal is calculated by filtering the input signal using a frequency 

27 domain filter with real coefficients and may be a time series or a set of spectral 

28 estimates. If the output is a time series then it may be passed to a digital to analog 

29 converter (DAC) and an analog anti-imaging filter to produce an analog output signal or 

30 it may be used as an input to subsequent signal processing. 

3! The estimator of the spectral components comprises four basic steps: 

32 1 . Calculation of the spectrum of the input signal. 
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1 2. Estimation of the signal and noise power in each frequency bin within the 

2 speech band (f_low->f_high Hz). 

3 3. Calculation of the gains (coefficients) of the frequency domain filter for each 

4 frequency bin, and 

5 4. Calculation of the spectral estimates by multiplying each input spectral 

6 component by the corresponding gain. 

7 This is basically the method of Graupe and Causey, and each of the processes is 

8 discussed below. 
9 

10 The estimates of the noise are updated during pauses in the information signal. 

11 These pauses are detected by looking at the power estimate to see if it exceeds a 

12 predetermined threshold, noise threshold, multiplied by noise[f] at each frequency. If 

13 the power estimate is above the calculated threshold then a thrsholdCnt[fJ is checked to 

14 see if it exceeds a predetermined value update_delay. 
15 

16 The spectral components of the input signal can be obtained by a variety of 

17 means, including band pass filtering and Fourier transformation. In one approach a 

18 discrete or fast Fourier transform is used to transform sequential blocks of N points of 

19 the input time series. A window function, such as a Harming window, can be applied, in 

20 which case an overlap of N/2 points can be used. A Discrete Fourier Transform (DFT) 

21 can be used at each frequency bin in the speech band or, alternatively, a Fast Fourier 

22 Transform (FFT) can be used over the whole frequency band. The spectrum is stored 

23 for each frequency bin within the speech band. For some applications it is desirable to 

24 have unequally spaced frequencies - in these applications a Fast Fourier transform 

25 cannot be used and each component may have to be calculated independently. In one 

26 approach the input spectrum, X, is calculated as the Fourier transform of the input time 

27 series, x, namely 

28 X= Fourier transform { x, window function, N }. 
29 

30 The power in the input spectrum is given by 

31 power - modulus squared {X}. 

32 
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Alternatively, a band pass filter may be used, in which case the power may be 
estimated by rectifying and smoothing the filter output. This version of a Graupe and 
Causey system is shown in Figure 1, Block Diagram 100. Input Time Signal 105 . x, is 
applied to a bank of band pass filters. One of these bandpass filters is represented by 
Bandpass Filter 110 in Figure 1. The output of Bandpass Filter 110 is Input Spectral 
Signal 115, referred to as X. The power of Input Spectral Signal 115 is measured by 
Input Spectral Signal Power Measurement 140, which generates Total Input Spectral 
Power Signal 165. The method requires that estimates be made for both Total Input 
Spectral Power Signal 165 and Noise Power Estimator Output 160. Noise Power 
Estimator Output 160 is generated by Noise Power Estimator 145 which utilizes a time 
constant related to the time over which the noise content of Total Input Spectral Power 
Signal 165 can be considered stationary. Total Input Spectral Power Signal 165 is 
estimated by Signal Power Estimator 155. From these estimates Wiener Gain 
Coefficients 170 is calculated by Wiener Gain Calculator 150, Wiener Gain Calculator 
150 determines the ratio of the power in the information signal, which is Total Input 
Spectral Power Signal 165, to the total power which is the sum of Noise Power 
Estimator Output 160 and Total Input Spectral Power Signal 165. For each frequency 
bin this is 

W = signal /( noise + signal ). 
In the method of Graupe and Causey the Wiener gain. W, is directly applied to the 
corresponding component of the input spectrum. In the unmodified scheme the spectral 
components of the output are given by multiplying Input Spectral Signal 115 by Wiener 
Gain Coefficients 170 in Multiplier 120. The result is 

r= W*X 

which is Output Spectral Signal 125. If Output Time Signal 135, y, is required it can be 
calculated by an inverse FFT (or DFT) and the 'overlap-add' method or by summing the 
components from individual channels using Channel Combiner 130. 

After each iteration k the output block of N time points is updated as 
yk(l:N) = inverse Fourier transform {Y,N} 
y k (l:N/2) = yk (l:N/2) + y k .j(N/2+l :N) 
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1 The first N/2 points of are then sent to Channel Combiner 130 or may be used for 

2 further processing. 
3 

4 An improved system is shown in Figure 2, Block Diagram 200. The additional 

5 features are described below. 
6 

7 Gain Modification 

8 Time Input Signal 205 is applied to Bandpass Filter 210. The output of 

9 Bandpass Filter 210 is applied to the input of Multiplier 220, and if a time signal output 

10 is desired Channel Combiner 230 is utilized to generate Time Output Signal 235. When 
ti the signal to noise ratio is low the direct use of the Wiener gain results in a residual 

12 noise which has a musical or artificial character. One improvement is the use of Gain 

13 Modifier 270, which reduces the musical nature of the residual noise. Gain Modifier 

14 270 receives inputs from Wiener Gain Calculation 250 and Noise Power Estimator 245. 

15 The output of Total Input Spectral Power Measurement 240 is also routed as an input to 

16 Gain Modifier 270. 
17 

18 Gain Modifier 270 is presented by Figure 3, Block Diagram 300. The 

19 instantaneous power of the information signal can be estimated as the product of the 

20 instantaneous power and the Wiener gain. This gives an estimate of the instantaneous 

21 signal to noise ratio, snr, in each frequency bin obtained by dividing Total Input 

22 Spectral Power 265 by Noise Power Estimator Output 260, which is accomplished by 

23 Divider 305, and using this quotient to modulate or multiply Wiener Gain Coefficients 

24 280. This is accomplished by Multiplier 325, and the output of Multiplier 325 is Signal- 

25 to-Noise Ratio Estimate 320. Hence 

26 snr — W * (power / noise). 

27 A function of the signal to noise ratio is then calculated by Function Modifier 315, and 

28 Modified Coefficients 275, which are denoted by the vector C, are calculated by 

29 dividing the output of Function Modifier 315 by the output of Divider 305. This is 

30 accomplished by Divider 310 and is done for each frequency, so that 

31 C = F{snr} * (noise /power) = F{snr) / (power / noise) 
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1 where F is a function of a single variable and is therefore well suited to implementation 

2 on a DSP as a look-up table or an analytic function. One form of the function F is 

3 given by 

{x ]n . x<snr0 
x + c, .v>w0- 

5 where c and snrO are constants. Other forms can used, but it is desirable that the 

6 ftmction is approximately linear at high signal to noise ratios. In particular the gain of 

7 Ephraim and Malah may be manipulated so that it can be implemented in this form. 
8 

9 Output Spectral Signal 225, Y, which is the estimate of the spectrum of the 

10 information signal, is calculated by multiplying 215 by the corresponding Modified 

1 1 Coefficients 275, as shown in Figure 2, so that for each frequency 

12 Y=C*X 
13 

14 Signal Estimation 

15 Ephraim and Malah describe a method for updating a signal to noise ratio. This 

16 method can be modified to give an estimate of Signal Power Estimator Output 285. 

17 Signal Power Estimator 255 uses the power in the output spectral signal Output Spectral 

18 Signal Power 290 which is calculated by Output Spectral Signal Power Measurement 

19 295 as shown in Figure 2. The method is shown in detail in Figure 4, Block Diagram 

20 400, and is given by 

21 sigl = maximum {power - noise, 0} 

22 sig2 = modulus squared {Y} 

23 signal = (I -beta) * sigl + beta * sigl 
24 

25 The difference between Total Input Spectral Power 265 and Noise Power 

26 Estimator Output 260 is calculated by Adder 405. The output of Adder 405 is half-wave 

27 rectified by Half-wave Rectifier 410. The output of Half-wave Rectifier 410 is Half- 

28 wave Rectifier Output Signal 415, and Half-wave Rectifier Output Signal 415 is 

29 weighted by (1-Beta) Weighting Function 420. Signal Power Estimator Output 285 is 

30 obtained as the sum of the output of (1-Beta) Weighting Function 420 and the output of 

31 (Beta) Weighting Function 430 by Adder 425. The output of (Beta) Weighting Function 
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1 . 430 is a weighted value of Output Spectral Signal Power 290. The weighting parameter 

2 beta used in the weighted sum is typically chosen to be greater than 0.9 and less than 1 . 
3 

4 Noise Estimation 

5 The estimates of the noise can be updated during the pauses in the information 

6 signal. The pauses can be detected by looking at the power estimate to see if it exceeds 

7 a predetermined threshold, noise threshold multiplied by noise[f] at each frequency. If 

8 the power estimate is above the calculated threshold then a thrsholdCnt[f] is checked to 

9 see if it exceeds a predetermined value updatedelay. If it does, the noise estimate is 

10 updated as 

1 1 Noise[f]=alpha*power+( 1 -alpha)*noise[f) 

12 Noise[f]=maximum {noise[f], minNoise} 

13 ThrsholdCnt[f]=0 

14 MinNoise is a constant that prevents noise[f] from being equal to zero. It is 

15 typically equal to 1 * 1 0exp-7. 

16 If a pause is not detected, thrsholdCnt[f] is incremented 

17 ThrsholdCnt[f]=thrsholdCnt[f]+ 1 
18 

19 This type of noise estimator is depicted in Figure 5, Block Diagram 500. Input 

20 Spectral Power 265 is applied to a first input of Adder* 575. Noise Power Estimator 

21 Output 260 is applied to the input of Time Delay 565. Time Delay 565 functions as a 

22 one-sample delay. The output of Time Delay 565 is multiplied by a constant, Noise 

23 Threshold, in Noise Threshold Multiplier 570. The output of Noise Threshold 

24 Multiplier 570 is routed to a second input of Adder 575. The output of Adder 575 is 

25 input to (>= 0?) Function 550. The inputs of the logical AND 556 enables Algorithmic 

26 Process 545 only if Function 550 and Function 555 are true. If the output of function 

27 556 is False Algorithmic Process 545 is disabled. The output of Logical AND 556 is 

28 applied to the input of Inverter 557. A false input at Inverter 557 will enable 

29 Algorithmic Process 540. 
30 

31 The sequence of First Algorithmic Process 540 will now be described. Input 

32 Spectral Power 265 is input to (alpha) Multiplier 505. The output of (alpha) Multiplier 
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1 505 is applied to a first input of Adder 510. The output of Adder 510 is an input to 

2 Multiplier 525. Time Delay 520 is a one sample delay. The output of Multiplier 525 is 

3 Noise Power Estimator Output 260. Noise Power Estimator Output 260 is applied to the 

4 input of Time Delay 565 and to the input of Time Delay 530. Time Delay 530 functions 

5 as a single sample delay. The output of Time Delay 530 is applied to the input of (1- 

6 alpha) Multiplier 515. The output of (1-alpha) Multiplier 515 is applied as a second 

7 input of Adder 510. 
8 

9 Second Algorithmic Process 545 produces an increment in the value of 

10 ThrsholdCnt, and is represented by (Increment thrsholdCnt) Function 535. 
ll 

12 Information Signal Detector 

13 T ne present invention operates to update estimates of the noise during pauses in 

14 the information signal. The presence of an information signal can be detected by 

15 looking at a weighted sum of the signal to noise components across frequency bins (a 

16 uniform weighting may be used). If this weighted sum is above a predetermined 

17 threshold, the signal is assumed to contain information and the noise estimate is 

18 updated. This is shown in Figure 6, Block Diagram 600. Signal-to-Noise Ratio 

19 Estimates 605 are weighted by Signal-to-Noise Ratio Weighting Coefficients 610 and 

20 then summed by Summer 615 to produce Summer Output Signal 630. S, before being 

21 input to Threshold Detector 620. The output of Threshold Detector 620 is Threshold 

22 Detector Output Signal 625. 

23 

24 One algorithmic example is described below: 

25 at each update number k 

26 X= Fourier transform { x, window Junction, N). 

27 FOR each frequency number fin speech band 
2 g power = modulus squared { X[f] } 

29 sigl = maximum {power - noise[f], 0} 

30 sig2 = modulus squared{ 17X7} 

31 signal = (I -beta) * sigl + beta * sig2 

32 W = signal /( noise [f] + signal ) 
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! snr = W * ( power /noise [fl ) 

2 C = F{snr} I ( power /noise [f] ) 

3 IF(power-noiseThreshold*noise[f]>=0 and thrsholdCnt<update_delay 

4 THEN thrsholdCnt[f]=thrshoIdCnt[f]+l 

5 OTHERWISE noise[fj=alpha*powert-(l -alpha)*noise[f] 

6 Noise[fj=max(noise[f], minNoise) 

7 ThrsholdCnt[f]=0 

8 ENDIF 

9 old _power[f] = power 

10 Y[f] = C*X[f] 

11 ENDFOR 

12 y/cCl 'N) = inverse Fourier transform {Y,N} 

13 ygl:N/2) - yk (l:N/2) + y k .j(N/2+l:N) 

14 At the end of each iteration, the signal ygl:N/2) provides an estimate of the 

15 information signal. If a pause is not detected in the information signal, then 

16 thrsholdCnt[f] is incremented: thrsholdCnt[fj = thrsholdCnt[f] + 1 . 
17 

18 As can be seen from the foregoing description the present invention teaches a 

19 method whereby a noise estimate may be calculated and utilized in an adaptive speech 

20 filter algorithm. The noise estimation method generates noise estimates only during 

21 pauses in the information signal, rather than continuously updating the noise estimates. 

22 This noise estimation and updating technique allows for faster convergence and quicker 

23 cancellation of interfering tones than prior art techniques. The algorithmic technique 

24 can be implemented on inexpensive digital signal processors. It typically will result in 

25 less processing time, and memory requirements are less. The method of the present 

26 invention avoids corruption of the noise estimates due to additive information signal 

27 content that is common in other methods of noise estimation. 
28 

29 While the invention has been particularly shown and described with reference to 

30 a preferred embodiment, it will be understood by those skilled in the art that various 

31 changes in form and detail may be made therein without departing from the spirit and 

32 scope of the invention. 
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l What is claimed is: 

2 

3 1. A method for estimating the frequency components of an information signal 

4 from an input signal containing both the information signal and noise, said method 

5 comprising: 
6 

7 producing a set of input frequency components, one for each frequency band and 

8 for each frequency component; 
9 

10 calculating the total power in each input frequency component; 

11 

12 estimating the power of the information signal included therein; 

13 

14 calculating a modified gain for each frequency band as a function of the total 

15 power, the estimate of the power of the information signal and an estimate of the noise 

16 power; 
17 

18 multiplying the input frequency component by said gain to thereby produce an 

19 estimate of the power of the frequency component of said information signal; 
20 

21 detecting a pause in the information signal, further comprising the following 

22 steps: 
23 

24 determining whether the estimate of the power of the frequency 

25 component of said information signal exceeds a first predetermined threshold, and 
26 

27 if the estimate of the power of the frequency component of said 

28 information signal exceeds the first predetermined threshold, then determining whether a 

29 threshold value exceeds a second predetermined threshold wherein if the pause is 

30 detected if the threshold value exceeds the second predetermined threshold, and 
31 
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1 updating the estimate of the noise power during the pause detected in the 

2 information signal. 
3 

4 2. A method as in claim 1 wherein if the threshold value does not exceed the 

5 second predetermined threshold, then incrementing the threshold value. 
6 

7 3. A method as in claim 1 and including the step of estimating the overall signal to 

8 noise ratio from a weighted sum of the estimated signal to noise ratios of each frequency 

9 band. 
10 

n 4. A method as in claim 3 and including the step of using said estimated overall 

12 signed to noise ratio to determine the presence of an information signal in the input 

13 signal. 
14 

15 5. A method as in claim 1 and including the step of recombining the estimates of 

16 the frequency components of said information signal to produce a noise reduced output 

17 signal. 
18 

19 6. A method as in claim 1 in which the power of the information signal is estimated 

20 from a combination of a previous estimate of the frequency components of said 

21 information signal and the positive difference between the power in the input frequency 

22 component and the estimate of the noise power. 
23 

24 7. A method as in claim 6, in which the gain in each frequency band is determined 

25 by the steps of: 
26 

27 estimating a Wiener gain from said estimate of the noise power and the estimate 

28 of the power of the information signed; 
29 

30 multiplying said Wiener gain by the ratio of the power of the input frequency 

31 component to the estimated noise power to produce an estimate of the signal to noise 

32 ratio; 
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1 

2 calculating a function of the estimated signal to noise ratio; and 

3 

4 dividing said function of the estimated signal to noise ratio by the ratio of the 

5 power of the input frequency component to the estimated noise power to thereby 

6 produce a modified gain. 
7 

8 8. A method as in claim 1 in which the filtering is performed via a Fourier 

9 transform. 
10 

n 9. a method as in claim 1 which is used as a preprocessor to a speech or voice 

12 recognition system. 

13 

14 10. A method as in claim 1 which is used for reducing noise in a communications 

15 system. 
16 

17 ii. The method of claim 1 wherein said means to produce a set of input frequency 

18 components comprises filtering the input signal through a set of band pass filters. 
19 

20 12. The method of claim 1 wherein said means to produce a set of input frequency 

2 1 components comprises calculating the Fourier Transform of the input signal . 
22 

2 3 13. A system for estimating the noise power of frequency components of an 

24 information signal from an input signal containing both the information signal and 

25 noise, said system comprising: 
26 

27 means to produce a set of input frequency components, one for each frequency 

28 band and for each frequency component; 
29 

30 a first calculating means for calculating the total power in each input frequency 

31 component; 
32 
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1 an estimating means for estimating the power of the information signal included 

2 therein and for updating an estimate of the noise power only during a pause detected in 

3 the information signal by the estimating means; 
4 

5 a second calculating means for calculating a modified gain for each frequency 

6 band as a function of the total power, the estimate of the power of the information signal 

7 and the estimate of the noise power; and 
8 

9 gain multiplying means for multiplying the input frequency component by said 

10 gain to thereby produce an estimate of the power of the frequency component of said 

1 1 information signal. 

12 

13 14. A system as in claim 13 wherein the estimating means comprises: 
14 

15 an adder that is provided with an input spectral power signal and a first 

16 predetermined threshold value; 
17 

18 a first comparison element that receives the estimate of the power of the 

19 information signal and the first predetermined threshold value from the first adder, 

20 wherein the first comparison element determines whether the estimate of the power of 

21 the information signal exceeds the first predetermined threshold value; and 
22 

23 a second comparison element coupled to the first comparison element that 

24 determines whether a threshold value exceeds a second predetermined threshold value if 

25 the estimate of the power of the information signal exceeds the first predetermined 

26 threshold value, wherein if the threshold value exceeds the second predetermined 

27 threshold value then the pause is detected and the estimating means updates the estimate 

28 of the noise power. 
29 

30 
31 
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1 15. A system as in claim 1 3 further comprising: 
2 

3 an increment element that increments the threshold value if the threshold value 

4 does not exceed the second predetermined threshold value. 
5 

6 16. A system as in claim 1 3 in which the second calculating means comprises: 
7 

8 means for estimating a Wiener gain from said estimate of the noise power and 

9 the estimate of the power of the information signal; 
10 

u Wiener multiplying means for multiplying said Wiener gain by the ratio of the 

12 power of the input frequency component to the estimate of the noise power to produce 

13 an estimate of the signal to noise ratio; 
14 

15 function calculating means for calculating a function of the estimated signal to 

16 noise ratio; and 
17 

18 division means for dividing said function of the estimated signal to noise ratio by 

19 the ratio of the power of the input frequency component to the estimate of the noise 

20 power to thereby produce a modified gain. 

21 

22 17. The system of claim 13 wherein said means to produce a set of input frequency 

23 components is capable of filtering the input signal through a set of band pass filters. 
24 

25 18. The system of claim 13 wherein said means to produce a set of input frequency 

26 components is capable of calculating the Fourier Transform of the input signal. 
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